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Preface 


This  thesis  is  the  third  in  a  series  of  theses 
which  have  studied  the  mechanism  by  which  organo-nitrogen 
compounds  present  in  fuels  of  origin  other  than  petroleum 
are  converted  to  NC>x  pollutants.  Due  to  current  uncertainty 
in  supply  of  petroleum  the  Air  Force  is  interested  in 
fuels  derived  from  coal  or  oil  shale.  A  concern  is  that 
use  of  these  fuels  carries  the  probability  of  high  levels 
of  NO  pollutants  will  be  formed  in  combustion. 

A 

The  pyrolysis  of  monomethylamine  was  studied  in  a 
single  pulse  shock  tube  facility  with  product  analysis  by 
gas  chromatography,  as  opposed  to  the  two  previous  theses 
which  used  in  situ  infrared  emission  techniques.  Results 
from  this  study  support  part  of  the  model  for  decomposition 
of  monomethylamine  developed  in  the  previous  two  theses, 
however,  further  work  needs  to  be  done.  An  activation 
energy  and  pre-exponential  factor  for  the  initial  step 
in  the  decomposition  reaction  was  determined. 

I  would  like  to  thank  those  people  whose  help  made 
completion  of  this  thesis  possible.  Especially,  I  would 
like  to  thank  Dr.  E.  A.  Dorko,  my  thesis  advisor  and 
instructor  at  the  Air  Force  Institute  of  Technology  (AFIT) 
who  provided  me  with  the  knowledge  and  guidance  to  complete 
this  thesis.  I  would  like  to  also  thank  Harold  L.  Cannon 
of  AFIT  for  his  help  in  operation  of  the  shock  tube,  and 
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William  Baker  for 'his  help  In  obtaining  needed  supplies. 
Finally,  I  want  to  thank  my  wife,  Lt.  Gail  M.  Nusz,  for 
-"-her  assistance  with  the  computer  work  and  plotting,  and 
for  her  typing  of  the  completed  thesis. 


William  R.  Nusz 
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Abstract 


The  decomposition  of  monomethyl amine  (MMA)  was  studied 
to  model  the  decomposition  mechanism.  The  decomposition  of 
a  2.5?  mixture  of  MMA  in  argon  was  studied  behind  a  reflected 
shock  by  use/ of  a  single  pulse  shock  tube  with  gas  chromato- 

,  rl  H 

graphy  of  the  products  CH^,  NH^,  and  unreacted  MMA..  The 
temperatures  behind  the  reflected  shock  ranged  from  727  °K 
to  1438  °K.  The  reaction  test  time  was  . 900  i/sec.  Total 
pressure  ranged  from  2.2  atm  to  8.03  atm.  The  rate 
constants  for  the  initiating  C-N  bond  scissure  were 

determined  by  a  parametric  study  using  a  numerical  model. 

M  r] 

The  concentrations  of  CH^ ,  NH^,  and  MMA  were  examined  to 
determine  the  reasonableness  of  the  model..  The  Arrhenius 


A 

re  is 


equation  for  the  initial  C-N  bond  scissure 'Is  lo^-^* 
exp(-23- 800/RT)  sec-1  as  determined  in  this  study. 
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A  SINGLE  PULSE  SHOCK  TUBE 


STUDY  OP  PYROLYSIS  OP  MONOMETHYLAMINE 

I.  Introduction 


Background 

The  Air  Force  is  presently  interested  in  the  production 
of  fuels  from  sources  of  origin  other  than  crude  petroleum, 
such  as  "synthetic"  crudes  produced  from  coal  or  from  oil 
shale.  One  problem  with  these  alternate  sources  is  that 
fuels,  specifically  JP-4  and  Jet  A,  produced  from  oil  shale 
or  coal  have  a  relatively  large  percentage,  .3  to  2%  by 
weight,  df  organo-nitrogen  compounds  (Ref  15:5),  as 
compared  to  less  than  .015?  for  petroleum  based  fuels. 

During  combustion  these  alternate  fuels  produce  large 
quantities  of  NOx  pollutants,  on  the  order  of  140%  more 
than  petroleum  based  fuels  (Ref  10:7).  The  symbol  NO 

*  X 

refers  to  the  three  oxides  of  nitrogen  NO,  N02 ,  and  N20^ 
which  are  usually  created  In  high  temperature  combustion 
(Ref  4).  N0x  compounds  in  the  atmosphere  participate  in 
the  formation  of  photochemical  smog  and  may  have  detrimental 
effects  on  the  stratospheric  ozone  layer  (Ref  4). 

The  Department  of  Transportation  has  recommended  basic 
combustion  research  to  reduce  the  amount  of  N0x  emission 
by ^aircraft.  Since  JP-4  and  Jet  A  produced  from  coal  and 
oil  shale  derived  crude  are  projected  to  be  produced  in  the 
future  (Ref  4),  the  Air  Force  is  interested  in  reducing 
the  amount  of  NO  formed  when  these  fuels  are  burned. 

X 


It  is  currently  perceived  that  in  combustion  modeling 
the  combustion  process  can  be  represented  by  finite  rate 
chemical  reactions  if  the  rate  data  are  known  accurately 
(Ref  7)-  Also  while  a  partial  mechanism  can  simulate  the 
heat  process  due  to  combustion,  knowledge  of  the  complete 
mechanism  is  necessary  in  order  to  predict  the  pollutants 
accurately  (Ref  7).  To  learn  how  to  reduct  the  amounts 
of  N0x  produced,  it  is  first  necessary  to' learn  the  mechanisms 
that  lead  to  formation  of  N0x  by  organo-nitrogen  compounds. 
Once  the  mechanism  is  understood,  the  critical  steps  that 
lead  to  N0x  formation  may  be  inhibited  by  adjustment  of 
reaction  conditions. 

For  the  case  of  combustion  of  petroleum-based  fuels, 
the  main  source  of  N0x  compound  is  known  to  be  fixation  of 
atmospheric  nitrogen.  The  reaction  is  thought  to  proceed  by 
the  Zeldovich  mechanism. 

Q2  +  M  -*•  20  +  M 

0  +  N2  -*■  NO  +  N 

N  +  02  +  NO  +  0 

It  is  known  that  this  mechanism  can  be  supressed  by  operating 
at  lower  temperatures.  However,  formation  of  N0x  from  organo- 
nitrogen  fuel  components  occurs  readily  even  at  low  combustion 
temperatures  with  a  conversion  rate  for  fuel-bound  nitrogen 


of  50-100$  (Ref  4).  This  knowledge  has  prompted  studies 
to  find  the  mechanisms  by  which  organo-nitrogen  compounds 
form  N0X.  A  compound  which  has  been  chosen  in  the  past 
as  a  suitable  model  for  such  studies  is  monomethylamine 
(MMA)  (Ref  12). 


Research  Objectives 

The  objective  of  this  thesis  was  to  make  a  kinetic 
study  of  the  decomposition  of  MMA  in  argon  as  a  step  toward 
modeling  the  mechanism  by  which  organo-nitrogen  compounds 
form  NOx-  The  pyrolysis  of  MMA  was  studied  in  a  single 
pulse  shock  tube.  The  products  of  decomposition  were  anal¬ 
yzed  by  vapor  phase  chromatography  (Ref  5).  The  arrhenius 
parameters  for  decomposition  of  monomethylamine  were  to  be 
calculated  by  matching  concentrations  of  CH^,  NH^,  and  MMA 
obtained  by  integration  of  a  numerical  model  to  the  con¬ 
centrations  of  the  same  species  found  experimentally. 


II.  Theory 


Mechanism 

Monomethylamine ,  CH3NH2,  was  chosen  for  this  study 
because  it  has  previously  been  used  for  kinetic  studies 
in  shock  tubes  using  infrared  emissions  techniques  (Ref  12, 
16).  MMA  is  the  simplest  organo-nitrogen  compound,  .is 

i 

stable,  readily  available,  and  is  analogous  to  ethane, 
which  has  a  well-known  decomposition. 

The  initial  steps  in  the  pyrolysis  of  MMA  has  een 
previously  studied.  The  most  recent  work  in  this  area 
gives • strong  evidence  to  support  the  following  mechanism 
(Ref  6) 


(1) 

CH3NH2  +  M 

CH3 

+ 

nh2  +  ] 

k„ 

(2) 

CH0  +  CH^NH- 
3  3  2 

CH^ 

+ 

CH-2NH2 

k-, 

(3) 

ch3  +  ch2nh2 

-.3 

ch4 

+ 

ch2=nh 

(4) 

nh2  +  ch3nh2 

k4 

-► 

NH3 

+ 

ch2nh2 

k 

(5) 

nh2  +  ch2nh2 

NH3 

+ 

ch2=nh 

(6) 

ch2=nh 

H2  ' 

f  : 

HCN 

This  mechanism  applies  well  for  times  from  0  to  ^350  ysec 
and  temperatures  from  1275  -  1400°K.  It  is  not  yet  known, 
however,  if  this  mechanism  is  general  enough  to  be  applic- 
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able  in  modified  form  in  a  lower  temperature  region  and  for 
longer  times.  At  lower  temperatures,  recombination  reactions 
are  likely  to  become  important.  To  account  for  recombination, 
the  following- reactions  should  be  added  to  the  above: 

C2H6 
CH3NH2 

A  study  by  Smith  &  Sawyer  (Ref  14)  using  a  flow  reactor 
for  the  pyrolysis  of  MMA  followed  by  analysis  of  the  products 
by  vapor  phase  chromatography  covered  a  temperature  range 
of  1087  -  1132  °K  with  a  residence  time  range  from  0  to  1.2x10 
sec.  They  reported  the  formation  of  N2,  solid  hexamethylene¬ 
tetramine  and  3-ethyl-4-methyl  pyrrole  in  addition  to  the 
products  of  the  mechanism  just  presented.  They  postulated 
that  the  hexamethylene-tetramine  and  pyrrole  were  formed 
from  the  CH2=NH  imine  intermediate  at  the  relatively  low 
temperatures,  while  at  higher  temperatures,  CH2=NH  forms 
H2  and  HCN . 

Kinetics 

Consider  the  simple  system: 

B  +  C 

D  +  E 

E  +  P 


(9) 

A 

kl 

-► 

k2 

(10) 

A  +  C 

-► 

k3 

(11) 

A  +  3 

(7) 

CH3  +  CH3 

(8) 

CH3  +  NH2 

CD 
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If  the  above  reactions  and  species  comprise  the  entire 
system  it  is  possible  to  write  the  chemical  mechanism  as 
.a  system  of 'linear  ordinary  differential  equations  in  terms 
of  the  species  concentration  and  the  rate  constants  as 
follows : 


(12) 

ft<A>  ■ 

-k1(A)  -  k2(A)(C)  - 

(13) 

5t<B>  ’ 

kx(A)  -  k3(A)(B) 

(14)  • 

3t(0)  " 

kx(A)  -  k2(A)(C) 

(15) 

tt<c>  ' 

k2(A)(C) 

(16) 

ft<E>  * 

k2(A)(C)  +  k3(A)(8) 

(17) 

!t<F>  - 

k3(A)(B) 

These  equations  may  then  be  evaluated  as  an  initial 
value  problem,  and  the  chemical  behavior  over  a  period  of 
time  may  be  found. 

The  eight  equations  of  the  mechanism  for  MMA  pyrolysis 
may  be  similarly  expressed  in  differential  equation  form  as 
follows : 

.  (18)  |t(CH3NH2)  =  -k1(CH3NH2)  -  k2(CH3NH2) (CH3) 

-k  (CH3NH2)(NH2)  +  kg(CH3)(NH2) 

(19)  |fc(CH2NH2)  =  k2(.CH3NH2)(CH3)  +  k^  (CH^H^  (NH2) 

-k3(CH2NH2)(CH3)  -  k4(CH2NH2)(NH2) 
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-k3(CH2NH2)(CH3)  -  k?(CH3)2 
-kg(CH3)(NH2) 


'  j  (21)  |t(NH2)  =  k1(CH3NH2)  -  k4(CH3NH2)(NH2) 

‘  ' 

i  |  -k5(CH2NH2)(NH2)  -  kg(CH3)(NH2) 

if  (22)  |t(NH3)  =  k4(CH3NH2)CNH2) 

l]  +k5(CH2NH2)(NH2) 

i  ’  * 

;  'i  (23)  ft(CH4)  =  Jc2(CH3NH2)(CH3) 

"  i 

+k3CCH2NH2)(CH3) 

kj  (24)  |t(CH2«NH)-  k3(CH2NH2)(CH3) 

•  +k5(CH2NH2)(NH2) 

|  (25)  ft(H2)  *  ft(HCN)  “  -k6(CH2=NH) 

|  Again,  this  system  may  be  treated  as  an  initial-value 

[  problem. 

:  The  rate  constants  for  these  reactions  are  temperature 

S  dependent  according  to  the  Arrhenius  equation 

(26)  k  -  Ae”Ea/RT 

!  In-  logarithmic  form  the  equation  is 

(27)  log  k  -  log  A  -  Ea/2.303RT 

l 

\  If  log  k  is  plotted  vs.  1/T,  a  straight  line  results  with 

I  a  slope  of  -Ea/2.303R  and  a  y-intercept  of  log  A. 


( 
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Shock  Tube  Studies 

Gas  phase  kinetics  have  been  studied  by  use  of  shock 
tubes  (Ref  8).  A  shock  tube  is  a  device  to  produce  a 
strong  planar  shock  wave  which  is  generated  by  bursting 
a  diaphragm  of  appropriate  material  which  separates  a 
section  of  high  pressure  gas  from  a  section  of  low  pressure 
gas.  The  driven  section,  at. low  pressure,  contains  the 
gas  to  be  reacted.  The  driver  section,  at  high  pressure, 
contains  an  inert  gas.  When  the  diaphragm  is  burst,  a 
shock  front  traverses  the  length  of  the  driven  section, 
both  heating  and  pressurizing  the  gas.  Upon  contacting 
the  end  wall  of  the  driven  section,  the  incident  shock 
is  reflected  back  into  the  already  heated  gas.  This  re¬ 
shocking  multiplies  again  the  pressure  and  temperature 
increase  previously  achieved.  Then,  in  a  single  pulse 
tube  as  used  in  this  study,  a  rarefaction  wave  caused  by 

gas  traveling  to  the  dump  tank  reaches  the  reaction  zone 

5  o 

and  cools  the  reacting  gas  at  about  5x10  K/sec,  effec¬ 
tively  "freezing"  the  chemistry  at  that  time. 

The  concentration  of  species  in  the  test  gas  can  be 
studied  in  several  ways.  The  method  of  interest  in  the 
current  study  involves  taking  a  sample  of  the  test  gas 
after  the  shock  and  analyzing  it  via  gas  chromatography 
(Ref  8). 

r  • 

Analysis  by  gas  chromatography  is  useful  because  it 
provides  information  on  the  absolute  quantities  of  end 
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products  of  the  reaction  under  study.  Gas  chromatography 

does  have  the  disadvantages  that  it  may  only  be  used  to 

determine  chemically  stable  species,  and  that  the  technique 

can  only  look  at  the  composition  after  one  particular 

experimental  time  period. 

The  quantities  required  to  determine  the  shock 

conditions  are  the  initial  pressure  of  the  driven  gas, 

P^,  the  initial  temperature  of  the  driven-  gas, 

the  composition  of  the  driven  gas ,  and  the  incident 

shock  velocity,  U.  (Ref  8).  The  reflected  shock  zone 

parameters  of  interest  are  pressure,  P._,  temperature, 

b 

Tp.,  and  the  residence  of  reaction  test  time,  t. which 

is  the  time  from  reflection  of  the  incident  shock  to  the 

arrival  of  the  rarefaction.  T^  and  P^  can  be  determined 

from  the  incident  shock  mach  number,  M,  where  M  =  U^/a^, 

a-^  =  speed  of  sound  in  the  gas,  and  y  is  the  heat  capacity 

ratio  =  C/C  .  The  ratios  of  Tr-A1..  and  Pr-/P,  are 
p  v  b  x  b  1 


h.  -  (2(y-1)M2 

T1 


(3-y))((3y-1)*t  -  2(v~l)) 
(y+1)2M2 


(29) 


il  m  (2yM2  -  ( Y-l) ? ( ( 3y-1)M2  -  2(y+1)) 
P1  (y+1)((y-1)M2  +  2) 


The  shock  parameters  were  calculated  using  the  NASA 
273  program  assuming  frozen  chemistry.  Program  inputs 
were  incident  shock  speed,  initial  temperature,  and 
thermodynamic  properties  of  the  test  gas  (Ref  13)-  The 
reaction  test  time,  t,  was  determined  experimentally. 
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Numerical  Integration 

The  differential  equations  of  a  chemical  system  are 
of  the  form: 

(30)  dy. 

w  =  fi  (yi  >vo>y?>  •  •  • 


dy, 

dt 


n' 


f,2(y1»y2*y3>-*  •  ,yn^ 


dy 

dt^  fn(y1>y2>y3> • • • >yn) 


subject  to  initial  conditions 


(3D  yx(o)  «  y1Q 


y2(°)  =  y 


20 


yn(0)  ■  yn0 

Now  consider  the  system  as  being  represented  by 
the  vector  equation: 

(32)  -  ?(y) 
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The  numerical  integration  formula  chosen  was 'the 


fourth  order  Runge-Kutta  type  which  is 


(33)  y1+1 


*1  +  “<K?<Si>  + 
*  S?<s|+i>> 


where 


(3*0 

yi 

+ 

ff(yi) 

(35) 

yi+h 

yi 

+ 

(36) 

yi+l  = 

yi 

+ 

*t?(y  Hh) 

is  evaluated  in  that  sequence  (Ref  11) .  By  input  of  the 
rate  constants  and  initial  species  concentration  for 
equations  18  to  23,  the  behavior  of  the  chemical  .system 
can  be  simulated.  The  numerical  model  should  behave 
as  the  experimental  model  given  the  same  conditions. 

By  comparing  the  results  of  the  numerical  integration 
against  the  experimental  results,  the  correct  rate  constants 
may  be  found. 


III.  Experimental 


Shock  Tube 

General  Description.  The  shock  tube  used  In  this  study 
was  constructed  of  stainless  steel  tubing.  The  tube  was  . 
of  1"  diameter  circular  cross  section  with  a  driver’ section 
54"  long  and  a  driven  section  72"  long.  A  dump  tank  of 
approximately  15  liter  capacity  was  attached  to  the  driven 
section  6.5  inches  downstream  of  the  diaphragm.  Figure  1 
shows  an  overall -layout  of  the  shock  tube.  The  vacuum 
system  consisted  of  an  Alcatel  Z2030C  roughing  pump  and 
a  Veeco  2-5W  oil  diffusion  pump  for  the  driven  section 
and  dump  tank.  With  this  system  the  driven  section  of 
the  shock  tube  could  be  evacuated  to  less  than  10  J  Torr. 


Dump  Tank 


Driver 


& 


K 


54"- 


To  Gas  Handling  and  Vacuum 


Driven  Section 


Diaphragm  location 
-  72"-- 


Figure  1,  Shock  Tube  Overall  Layout 


Test  Gas.  The  test  gas  mixture  of  2.5#  MMA  was 
prepared  from  MMA  (98.0#  pure)  and  argon  (99.995#  pure), 
both  purchased  from  Matheson  Gas  Products  and  used  with 
no  further  processing.  The  sample  was  prepared  in  a  stain¬ 
less  steel  bottle  by  the  following  procedure:  First, 
the  bottle  was  evacuated  to  a  pressure  of  less  than  10-^ 
Torr,  then  filled  with  MMA  to  a  pressure  that  would  give 
the  desired  concentration  of  MMA  at  the  final  total 
pressure.  The  bottle  was  then  closed  and  the  gas  manifold 
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pumped  down  to  less  than  10  Torr.  The  manifold  was  then 
pressurized  with  argon  to  at  least  double  the  pressure  of 
MMA  in  the  bottle  to  prevent  excape  of  the  MMA,  then  the 
valve  to  the  bottle  was  opened  and  the  bottle  filled  with 
argon  to  the  desired  total  pressure.  The  gas. sample  was 
allowed  to  stand  overnight  before  use  to  insure  homogeneity. 

Instrumentation .  The  purpose  of  the  instrumentation 
on  the  shock  tube  was  to  provide  accurate  information  on 
initial  driven  section  pressure,  incident  shock  velocity, 
and  dwell  time  of  the  reflected  shock  at  the  back  wall  of 
the  tube.  Pressures  from  1  to  200  Torr  were  measured  by 

a  Wallace  &  Tiernan  FA-l6o  aneroid  type  pressure  gauge. 
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Pressures  in  the  range  of  1  to  10  Torr  were  measured 
by  a  Consolidated  Vacuum  Corporation  DV-180  cold  cathode 
gauge.  Incident  shock  velocity  was  measured  by  two 
Kistler  piezo-electric  pressure  transducers,  placed  10 


inches  apart  as  shown  in  Figure  2.  These  transducers  were 
connected  through  charge  amplifiers  to  an  electronic  counter. 
The  time  interval  was  measured  by  a  Hewlett-Packard  5325B 
counter  with  'accuracy  of  ^0.1  ysec.  The  reaction  rest 
time  of  the  reflected  shock  was  measured  by  a  Kistler 
piezo-electric  transducer!  located  in  the  back  wall  of  the 
tube,  and  connected  to  a  Hewlett-Packard  551  oscilloscope 
through  a  charge  amplifier.  The  oscilloscope  was  triggered 
by  the  same  transducer  which  started  the  timer  to  measure 
incident  shock  velocity.  Figure  3  shows  a  sample  oscilloscope 
trace  reflected  shock  dwell-  time. 

Shock  Tube  Operation.  The  shock  tube  was  operated 
by  first  evacuating  all  sections  of  the  tube.  The  dump 
tank  was  then  filled  with  argon  to  the  desired  initial 
pressure  and  isolated  from  the  driven  section  which  was 
then  again  evacuated.  The  driven  section  was  then  filled 
with  the  test  gas  mixture  and  shocked  by  bursting  the 
diaphragm  by  admitting  high  pressure  helium  to  the  driven 
section.  Diaphragm  material  was  1.6  mil  mylar  sheet 
weakened  by  scoring  with  a  divider  point  to  allow  bursting 
in  the  proper  pressure  range.  A  more  detailed  description 
of  the  operation  appears  in  Appendix  A. 

Vapor  Phase  Chromatography 


V  •  Description.  Analysis  of  the  product  gas  mixture  was 
performed  on  a  Hewlett-Packard  model  720  gas  chromatograph 
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equipped  with  thermal  conductivity  detector.  Column 
materials  initially  considered  were  Chromosorb  102, 
Chromosorb  10 3 ,  Chromosorb  104,  Porapak  P,  Porapak  N, 

Linde  Molecular  Sieve  13-X,  and  Linde  Molecular  Sieve  5A. 
Chromosorb  103  was  found  to  be  suitable  for  use  in  deter¬ 
mining  MMA  and  NH^  concentrations,  and -Linde  Molecular 
Sieve  5A  was  found  to  be  suitable  for  the  determination 
of  CHjj.  The  operating  conditions  of  the  columns  are 
listed  in  Table  I.  Typical  chromatograph  traces  for  each 
column  are  reproduced  in  Figure  4.  The  Sargent  Model 
SR  x-y  recorded  used  for  recording  the  traces  from  the 
chromatograph  was  equipped  with  a  disc  integrator  which 
was  used  to  determine  the  areas  under  the  peaks.  The 
chromatograph  was  calibrated  using  a  mixture  of  5%  MMA, 

5%  NH^,  and  5%  CH^  in  argon.  The  calibration  curves 
appear  in  Figures  5-7. 

Operation.  At  the  start  of  each  day  both  columns 
were  seasoned  by  the  injection  of  a  2.0  ml  sample  of  the 
5%  MMA,  5%  NH^,  and  5%  CH^  mixture  at  ^350  Torr.  Samples 
were  introduced  by  a  Precjsion  Sampling  2 . 0  ml  gas-tight 
syringe  and  the  pressure  of  the  sample  in  the  gas  bottle 
was  determined  by  use  of  a  mercury  manometer.  The 
attenuation  of  the  chromatograph  amplifier  was  increased 
to  allow  the  argon  peak  to  remain  within  the  range  of  the 
x-y  recorder  and  reset  to  maximum  sensitivity  after  passage 
of  the  argon  peak.  The  sample  size  and  pressure  were 
recorded  on  the  chromatograph  trace. 


Table  I 


Column  Operating  Conditions 

Column:  V*  O.D.  Stainless  Steel,  6' 
Packing:  Chromosorb  103  80/100  mesh 

Temperature  Cycle:  150°C,  Isothermal 
Carrier  Gas  Flow:  40  std  ml/min  He 
Carrier  Gas  Pressure:  35  PSIG 
To  Detect:  CH3NH2>  NH^ 


Column:  V*  O.D.  Stainless  Steel,  6* 
Packing:  Linde  Molecular  Sieve  5A 
Temperature  Cycle:  50°C,  Isothermal 
Carrier  Gas  Flow:  40  std  ral/min  He 
Carrier  Gas  Pressure:  35  PSIG 
To  Detect:  CH^ 


VPC  CALIBRATION  PLOT  FOR  NH3 
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FIGURE  6,  NH3  CALIBRATION  CURVE 
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Data  Reduction 

Chromatograph .  For  each  sample  Introduced  Into  the 
chromatograph,  the  initial  number  of  moles  of  MMA  present 
bef.ore  the  shock  was  calculated  from  measurement  of  the 
pressure  and  volume  of  the  sample  and  application  of  the 
ideal  gas  equation  of  state.  The  area  under  each  of  the 
species  peaks  on  the  chromatograph  traces  was  determined 
by  the  disc  integrator  and  the  number  of  moles  of  each 
species  (MMA,  CHjj ,  and  NH^)  was  then  determined  from  the 
calibration  curves.  The  concentration  of  each  species 
present  was  then  expressed  in  the  form  of  (moles  of 
species )/( 100  initial  moles  of  MMA). 

Shock  Tube.  The  initial  composition,  pressure, 
temperature,  and  the  incident  shock  velocity  were  used 
as  inputs  for  the  NASA  273  program.  The  shock  parameters 
were  computed  using  the  assumption  of  frozen  chemistry. 

Over  the  temperature  range  in  this  study,  the  reaction 
test  time  is  dependent  primarily  upon  tube  geometry  and 
is  900  usee. 

Numerical  Smoothing.  The  final  species  concentrations 
were  plotted  as  a  function  of  temperatures  and  subjected  to 
a  fourth  order. least  squares  polynomial  curve  fit  using 
the  CC6600  library  subprogram  PLSCF  (Ref  1).  The  program 
is ^listed  in  Appendix  B.  The  least  squares  polynomial 
curves  for  decomposition  of  MMA  and  formation  of  CH^  and 


NH^  thus  obtained  were  evaluated  at  50  K  intervals  from 
1000°K  to  1450°K.  These  evaluations  were  then  taken  to 
be  the  final-  species  concentration  for  the  given  temperature 
for  further  processing. 

Numerical  Integration.  Equations  18  to  24  were 
incorporated  into  the  numerical  integration  program  which 
is  listed  in  Appendix  C.  The  program  used  the  CC660-0 
library  integration  subprogram  RKDP  (Ref  1).  Inputs  to 
the  program  were  initial  concentrations  of  each  species 
and  all  rate  constants  except  the  particular  rate  constant 
that  was  being  determined.  Initial  assumed  rate  constants 
appear  in  Table  II.  The  program  worked  by  assuming  an 
upper  and  lower  limit  for  the  rate  constant  under  search, 
calculating  the  species  concentrations  in  each  case, 
comparing  the  numerical  results  with  the  experimental, 
then  halving  the  interval  and  repeating  the  process  until 
the  rate  constant  is  known  to  the  desired  tolerance. 

was  determined  from  1000°K  to  1400°K.  The  k's 
were  then  plotted  as  log  k  vs.  1/T  and  the  appropriate 
Arrhenius  parameters  were  determined.  The  time  interval 
over  which  the  system  was  integrated  was  chosen  as  that  time 
between  the  arrival  of  the  reflected  shock  and  the  time 
that  pressure  behind  the  reflected  shock  had  fallen  to 
3/4  of  the  maximum  pressure  attained. 
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•  Table  II 

'  Initially  Assumed  Arrhenius  Parameters 


k 

log  A 

Ea(kcal/mole) 

Source 

kl 

17. 

79. 

(Ref  3:70) 

^2  *  ^4 

11. 

8.7 

(Ref  9) 

k3,k5 

10. 

10.0 

(Ref  9) 

k7 

10.5 

0. 

(Ref  3:104) 

kfl 

10.5 

0. 

(Ref  3:104) 
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IV.  Results 

For  each  of  the  end  product  species,  MMA,  CH^,  and 
NH^  that  were  analyzed,  plots  of  amount  of  species  present 
after  shock  vs.  shock  temperature  were  plotted  and  a  least 
squares  polynomial  curve  fit  as  described  earlier  was  made. 

For  MMA,  there  was  essentially  no  decomposition  at 
shock  temperatures  below  1000°K,  and  an  evenly  increasing 
rate  of  decomposition  with  increasing  temperature  until  at 
shock  temperatures  greater  than  1400°K  there  was  no  detec¬ 
table  amount  of  MMA  in  the  gas  sampled  after  shock.  The 
plot  for  MMA  appears  in  Figure  8. 

Essentially  no  NH^  was  formed  at  shock  temperatures 
below  1000°K,  with  an  increasing  yield  to  about  1275°K 
followed  by  a  rapid  fall-off  with  essentially  no  remaining 
NH^  for  shock  temperatures  above  1400°K.  The  rise  of  NH^ 
concentration  corresponds  directly  to  the  decomposition 
of  MMA,  the  fall-off  at  higher  temperatures  is  not  readily 
explainable  by  the  chosen  mechanism.  The  plot  for  NH^ 
appears  in  Figure  9. 

The  formation  of  CH^  was  also  .at  variance  with  the 
postulated  mechanism.  Essentially  no  CH^  was  detected 
for  temperatures  of  less  than  1200°K,  then  a  strongly 
increasing  yield  of  CH^  was  noted  for  temperatures  up  to 
the* high  limit  of  1400°K.  The  plot  for  CH^  appears  in 
Figure  10. 
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The  rate  constants  for  that  were  found  by  program 
KFIND  are  plotted  as  log  k  vs.  1/T  in  Figure  11.  The 
activation  energy  and  pre-exponential  factor  found  by 
this  plot  differ  greatly  from  expected  values. 

The  concentration  of  MMA,  NH^,  and  CH^  as  determined 
by  least  squares  curve  fit  are  compared  with  the  concentrations 
of  those  species  as  predicted  in  the  numerical  model  in 
Figures  12,  13,  and  14.  Note  that  the  calculated  concentration 
from  the  numerical  model  for  MMA  closely  follows  the 
observed,  the  calculated  concentration  vs.  the  observed 
concentration  for  NH^  is  in  good  agreement  to  1250°K,  and 
the  calculated  concentration  of  CH^  is  far  greater  than 
that  observed,  and  indeed,  methane  formation  is  not  noted 
at  temperatures  below  120Q°K. 


26 


600. 


700.  800.  900. 

TEHP.  DEGREES  K 


2000.  2200.  2200.  2300 


FIGURE  9,  CONCENTRATION  DATA  AND  LEAST  SQUARES  FIT  FOR  NH3 
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FIGURE  U.  PLOT  OF  LOG  K1  VS.  1/T  TO  DETERMINE  LOG  R 
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FIGURE  12,  CAPARISON  OF  EXPERIMENTAL  AND  CALCULATED  MMA  CONCENTRATION 
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FIGURE  13.  COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED  W3  CONCENTRATION 


HOLES/ 100  JNIT 


/ 


TEMP.  DECREES  K 


FIGURE  14,  COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED  CH4  CONCENTRATION 
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V.  Discussion 

The  most  unusual  result  obtained  in  this  study  is 
the  low  values  for  the  Arrhenius  parameters  obtained 
for  the  C-N  bond  scissure  of  MMA.  The  Arrhenius  para¬ 
meters  are  compared  in  Table  III.  Values  of  k^  lower 
than  those  that  would  be  expected  theoretically  have  been 
reported  elsewhere  (Ref  6,14)  and  have  been  attributed  to 
a  free  radical  mechanism,  such  as  the  postulated  mech¬ 
anism  used  here. 


Table  III 

Arrhenius  parameters  for  MMA  Decomposition 


log  A 

Ea(kcal/mole) 

Source 

14.62 

43-3 

(Ref  12) 

10.84 

48.15 

(Ref  5) 

17.0 

79.0 

(Ref  3) 

9.0 

58.6 

(Ref  14) 

6.64 

23.8 

this  work 
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formation  of  GH^  lagged  far  behind  that  expected.  Smith 
&  Sawyer  (Ref  14)  also  observed  a  rate  of  formation  of 
.CHjj  only  .42  times  the  rate  of  formation  of  NH^. 

As  an  experiment  in  this  study,  the  Arrhenius  para¬ 
meters  for  kg  and  k^  were  changed  by  decreasing  the  pre¬ 
exponential  factor  by  two  orders  of  magnitude  and  increasing 
Ea  by  50%  in  an  effort  to  supress  CH^  formation  in  t’he  numerical 
model.  This  change  did  not  supress  CH^  formation,  which 
is  the  expected  result,  as  the  reactions  involving  the 
radical  species  are  rapid  enough  to  allow  for  a  quasi¬ 
equilibrium  situation  with  very  low  CH^  and  NHg  concen¬ 
trations;  It  is  also  of  interest  that  due  to  the  low 
radical  concentration,  the  recombination  reactions  producing 
CH^NHg  and  CgHg  are  not  significant  as  these  reactions  are 
second  order  in  the  concentrations  of  the  radical  species. 

The  flow  reactor  study  of  Smith  &  Sawyer  (Ref  14) 
reported  other  results  which  are  interesting  to  compare 
to  the  results  obtained  here.  At  the  1060  to  1140°K 
temperature  range  they  operated  at,  large  quantities  of 
residue  collected  on  the  sampling  probe  used.  The  two 
main  constituents  of  the  residue  were  hexamethylenetetramine 
and  3-ethy 1-4-methyl  pyrrole.  The  pyrrole  could  be  formed 
by  an  attack  of  the  CHg=NH  species  by  the  CH^  species  through 
a  series  of  reactions  and  could  account  for  the  relative 
absence  of  CH^  at  temperatures  below  1200°K.  Above  1200°K, 


the  CH2=NH  tends  to  shed  hydrogen  and  become  H2  and 
HCN  as  in  equation  6. 

The  disappearance  of  NH^  at  temperatures  above 
1300°K  is  somewhat  of  a  puzzle.  The  NH^  is  probably 
not  dissociating  and  recombining  as  N2  and  H2,  as 
Bamford  (Ref  2:12)  reports  that  no  significant  homogeneous 
dissociation  of  NH^  occurs  below  about  2000°K.  The  same 
source,  however,  reports  that  heterogeneous  (i.e.  at 
vessel  wall)  decomposition  may  occur  at  lower  temperatures. 
It  can  be  assumed  that  NH^  is  reacting  with  other  species 
in  the  gas  mixture  at  temperatures  above  1275°K,  rather 
than  undergoing  homogeneous  decomposition. 
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VI .  Conclusions  and  Recommendations 

Conclusions 

As  a  result  of  this  study,  the  following  conclusions 
are  made: 

1.  The  postulated  mechanism  of  a  C-N  bond 
sclssure  followed  by  a  series  of  free  radical  reactions  is 
a  good  description  for  the  process  of  decomposition  of 
MMA  and  formation  of  NH^. 

2.  The  postulated  mechanism  is  not  a  good 
model  for  formation  of  CH^  or  for  tracking  of  the  carbon 
containing  species. 

3-  As  determined  by  the  methods  of  this  study, 
the  Arrhenius  equation  for  C-N  bond  scissure  of  mono- 
methylamlne  is 

k  «  106-6W-23-800/RT)  sec’1 

4.  The  numerical  model  behaves  in  a  consistent 
fashion  over  the  range  of  temperatures  and  concentrations 
in  this  study. 

Recommendations 

To  help  better  answer  some  of  the  problems  suggested, 
the  following  recommendations  are  made : 

1.  Procedures  for  better  determination  of  gas 
concentration  by  gas  chromatography  be  developed,  with 
special  attention  being  given  to  resolution  of  Ng, 

NgHjj,  C2Hg,  C2H^,  C2H2,  HCN,  and  Hg.  Use  of  a  flame 
ionization  detector  for  the  species  containing  carbon 
would  greatly  increase  sensitivity.  Other  analytical 
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techniques  such  as  mass  spectrometry  might  also  be  employed. 

2.  Theoretical  study  of  Cf^NH  and.  CH^  reactions 
should  be  made  to  attempt  to  model  the  behavior  of  these 
species  for  incorporation  into  the  numerical  model. 

3.  Further  numerical  studies  should  be  made 
incorporating  the  reactions  developed  in  recommendation 
2  to  try  to  better  approximate  the  experimental  behavior 
Of ;the  monomethylamine  pyrolysis  thus  developing  a  better 
understanding  of  the  decomposition  mechanism. 


ASD  Computer  Center,  "Subprogram  Library  Guide". 
Revision  E,  September,  1978. 

Bamford,  C.  H.,  C.  F.  H.  Tipper.  Comprehensive 
Chemical  Kinetics .  Volume  4.  Amsterdam: 
Elsevier  Publishing  Co.,  1972. 

Benson,  Sidney  W.  Thermochemical  Kinetics.  New  York: 
John  Wiley  &  Sons ,  Inc . ,  1968 . 

Blazowski,  William  S.  "Aircraft  Altitude  Emissions: 
Fundamental  Concepts  and  Future  R  &  D 
Requirements".  AIAA  Paper  75-1017,  AIAA  1975 
Aircraft  System  and  Technology  Meeting,  Los 
Angeles,  California,  August  4-7,  1975. 


Burcat,  A.,  G.  Skinner,  R.  Crossley  and  K.  Scheller. 

"High  Temperature  Decomposition  of  Ethane". 
International  Journal  of  Chemical  Kinetics,  V, 
345-352  (1973). 

Dorko,  E.  A.,  et^  al.  "Initial  Shock  Tube  Studies  of 

Monomethylamine" .  Journal  of  Physical  Chemistry, 

81 :  297  (1979). 

Garvin,  D.,  R.  L.  Brown,  R.  F.  Hampson,  M.  J.  Kuylo 
and  W.  Tsang.  "Summary  Report  on  the  Workshop 
on  High  Temperature  Chemical  Kinetics: 

Applications  to  Combustion  Research".  National 
Measurement  Laboratory,  National  Bureau  of 
Standards,  Washington  D.C.,  December  1978. 

Glass,  I.  I.  Shock  Tubes.  UTIA  Review  No.  12,  Part  II 
Institute  of  Aerophysics,  University  of  Toronto, 

May,  1958. 

Gray,  Peter  and  J.  C.  J.  Thynne'.  "Hydrogen  and  Deuterium 
Atom  Abstraction  from  Methylamine  and  Deuterated 
Methylamines" .  Trans.  Faraday  Soc.,  59  ••  2275 
(1963). 

Hardin,  M.  C.  "The  Combustion  of  Shale  Derived  Marine 

Diesel  Fuel  at  Marine  Gas  Turbine  Engine  Conditions" 
1976  Spring  Technical  Meeting  of  the  Central 
States  Section  of  the  Combustion  Institute, 

Columbus,  Ohio,  April  5-6,  1976. 

Hornbeck,  Robert  W.  Numerical  Methods.  New  York: 

Quantum  Publishing.  1975. 


Pchelkin,  Nicholas  R.  Initial  Shock  Tube  Kinetic 
Studies  of  Monomethylamlne .  GEP/CH/76D-1 ; 
Dayton:  Air  Force  Institute  of  Technology, 
December  1976. 


13.  Robinson,  P.  T.  and  K.  A.  Holbrook.  Unimolecular 

Reactions.  London: .  Wiley-Interscience,  1972. 

14.  Smith,  0.  I.  and  R.  F.  Sawyer.  The  Thermal  Pyrolysis 

of  Methylamine.  UCB-ME-76-2;  Berkeley:  University 
of  California;  April,  1976.  • 

15.  Vogt,  R.  A.  and  N.  M.  Laurendea'u.  ”N0x  Formation  from 

Coal  Nitrogen:  A  Review  and  a  Model”.  1976 
Spring  Technical  Meeting  of  the  Central  States 
Section  of  the  Combustion  Institute,  Columbus, 
Ohio,  April  5-6,  1976. 

16.  Wert,  John  C.  Further  Shock  Tube  Kinetic  Studies  of 

Monomethylamlne .  GEP/AA/77D-1;  Dayton :  Air 
Force  Institute  of  Technology,  December  1977. 


40 


/ 


APPENDIX  A 

SHOCK  TUBE  OPERATING 
PROCEDURE 


n 


The  shock  tube  was  operated  by  use  of  the  following 
procedure : 

1.  A  diaphragm  and  sample  bottle  were  installed 
and  all  sections  of  the  tube  were  evacuated  with  the  driven 
section  and  dump  tank  evacuated  to  less  than  10-^  Torr. 

2.  The  driven  section  and  dump  tank  were  filled 
with  argon  to  the  desired  initial  pressure,  then  the  dump 
tank  was  isolated  from  the  driven  section; 

3-  The  driven  section  and  sample  bottle  were 
evacuated  to  10  J  Torr  or  less  and  the  sample  bottle  was 
then  isolated  from  the  system. 

4.  The  driven  section  was  filled  with  test 

gas  to  a  pressure  2-3  Torr  higher  than  the  desired  initial 
pressure.  This  was  to  ensure  that  the  pure  argon  in  the 
dump  tank  would  not  flow  into  the  driven  section  and 
change  the  test  gas  mixture  when  the  driven  section  and  dump 
tank  were  connected. 

5.  The  counter  and  oscilloscope  were  reset, 
and  the  shutter  on  the  oscilloscope  camera  was  opened. 

6.  The  dump  tank  was  reconnected  to  the  driven 
section  by  opening  the  valve  between  them.  The  driven 
section  was  then  isolated  from  the  gas  manifold  to  prevent 
damage  to  the  pressure  gauges . 

7.  High  pressure  helium  was  introduced  to  the 
driven  section  until  the  diaphragm  was  burst  by  the  pressure. 
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8.  A  sample  of  the  gas  after  the  shock  was 
taken  by  opening  the  previously  evacuated  sample  bottle 
to  the  driven  section  of  the  tube  for  10  seconds. 

9.  The  counter  reading  was  recorded  and  the 
photograph  of  the  oscilloscope  developed. 

10.  The  tube  was  pressurized-  with  argon  to 
atmospheric  pressure  before  opening  to  change  diaphragms. 


APPENDIX  B 

LEAST  SQUARES  CURVE  FIT 
PROGRAM 
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PROGRAM  CFIT( INPUT, OUTPUT .PLOT) 

100  F0RMAT(F4.0,8X,F4.1) 

200  •  FORMAT ( "  "  ,IER=" ,12  ,  "NDEG=" ,12  ,  "NMAX=" ,12) 

300  FORMAT ( "  ",C(" ,11,")  ",F12.6) 

400  FORMAT ( "  ",F6.0,2X,F6.1) 

DIMENSION  X(49),Y(49),W(49),C(7) , WORK ( 200 ) 

DIMENSION  XP( 51) ,YP(51) ,U(22) ,V(22) 

DO  10  1=1,49 

10  READ  100, X(I) ,Y(I) 

DO  15  1=1,49 

15  PRINT  400,X(I),Y(I) 

NP=49 

NDEG=4 

NMAX=0 

IN=0 

W(1)=0.0 

CALL  PLS CF ( X , Y , W , NP , NDEG , NMAX , C , IN , XD , XO , WORK , IER ) 

PRINT  200, EIR, NDEG, NMAX 
NDEG=NDEG+1 
DO  20  1=1, NDEG 

20  PRINT  300,I,C(I) 

DO  30  1=1,20 

U(I)=950.  +  (I*25.  ) 

V(I)=(U(I)**4)»C(5)+(U(I)**3)*C(4)+(U(I)**2)C(3)+ 

1  (U(I)*C(2))+C(1) 

30  PRINT  400,U(I) ,V(I) 

DO  40  1=1,49 

xp(D-xd) 

40  YP ( I ) =Y( I ) 

CALL  PLOT(l. ,1. ,-3) 

CALL  SCALE(XP,6. ,49,1) 

CALL  SCALE(YP,6. ,49,1) 

U(21)=XP( 50) 

U(22)=XP(51) 

V(21)=YP(50)  • 

V(22)=YP(51) 

CALL  AXIS ( 0 . , 0 . , 15HTEMP .  DEGREES  K.-15.6 . ,0 . ,XP(50) ,XP(51) ) 
CALL  AXIS ( 0 . ,0 . , 24HMOLES/100  INIT. MOLES  MMA,24 ,6 . ,90, 

1  YP(50) ,YP(51) ) 

CALL  LINE(XP, YP, 49, 1,-1, 2)  . 

CALL  LINE(U,V, 20, 1,0,0) 

CALL  PLOTE(N) 

STOP 

END 


PROGRAM  KPIND( INPUT, OUTPUT) 

100  FORMAT(F4.0,lX,F7.5,2X,F4.1) 

200  FORMAT ( "  ", " IER= ", 12 , 1 X , "TEST® " , F6 . 2 , 2 X , " COUNT® " , F4 . 1 ) 

300  FORMAT ( "  ","TEMP=",F6 .1,2X,"K=",F13.1) 

EXTERNAL  F 
COMMON  RK(8) 

DIMENSION  Y(8),P(6),TG(9) , CONI (9 ) ,CON(9 ) 

DO  10  1=1,9 

10  READ  100 ,TG(I) ,CONI(I) ,CON(I) 

DO  20  NTRY=1,9 
C  COMPUTE  RATE  CONSTANTS 

RK(l)=10.*»17.*EXP(-79000./(r.98?«TG(NTRY))) 

RK(2)=10. **ll.*EXP(-8700./(1.98?»TG(NTRY))) 

RK(4)=RK(2) 

RK(3)®10.**11.*EXP(-10000./(1.987*TG(NTRY))) 

RK(5 )=RK( 3) 

RK(7 )®10 • **10 . 5 
RK(8)=RK(7) 

CT=0 . 

RU=RK(1)*8. 

RL=RK(1 ) * . 2 

1  RM= (RU+RL) /2 . 

CT=CT+1 . 

RK(1)=RM 

X=0 . 

DX= .0009 
N=6 

DO  2  1=1,6 

2  Y(I)=0 . 

Y ( 1 ) ®  CONI ( NTRY ) 

CALL  RKDF(X,Y,N,DX,IER) 

TEST  (Y(l)  CONI (NTRY)) *100. 

PRINT  200, IER, TEST, CT 
IF ( TEST. LT. CON (NTRY) )RU=RM 
IF ( TEST . GT . C  ON ( NTRY ) ) RL=  RM 
IF ( CT . GT . 2 5 • ) GO  TO  20 

IF(ABS(  (TEST-CON (NTRY)  )/?CON(NTRY) )  .GT.  ,01)GO  TO  1 
20  PRINT  300 ,TG(NTRY) ,RK(1) 

STOP 

END 


47 


SUBROUTINE  F(X,Y,P,N) 

COMMON  RK( 8) 

DIMENSION  Y(N) >P(N) 

P(1)--RK(1)*Y(1)-RK(2)*Y(1)*Y( 3)-RK(4)«Y(l)#Y(4) 

1  +RK(8)*Y(3)*Y(4) 

P(2)»RK(2)*Y(1)*Y(3)+RK(4)*Y(1)*Y(4)-RK(3)*Y(2) 

1  * Y ( 3 ) -RK (5)*Y(2)#Y(4) 

P(3)-RK(1)*Y(1)-RK(2)*Y(1)*Y(3)-RK(3)*Y(2)«Y(3) 

1  -RK( 7 ) *Y ( 3 ) *Y ( 3 ) -RK ( 8 ) *Y ( 3 ) *Y ( 4 ) 

P(4)»RK(1)«Y(1)-RK(4)*Y(1)*Y(4)  -  RK(5)*Y(2)*Y(4) 

1  RK(8)*Y(3)*Y(4) 

P  (  5 )  *RK  (  2.)  * Y  ( 1 )  *Y  ( 3  )  +RK  ( 3 )  *Y  ( 2 )  * Y  (  3 ) 
P(6)=RK(4)*Y(1)*Y(4)+RK(5)*Y(2)*Y(4) 

RETURN 

END 
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APPENDIX  D 


EXPERIMENTAL  DATA 


1153  .  70.6 
1111  91.8 
1054  93.7 
1190  40.7 
1139  57.6 

84.4 
59.2 
0.0 
0.0 
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